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We develop a new approach to study the nonlinear evolution in the large-scale structure of the Universe both 
in real space and in redshift space, extending the standard perturbation theory of gravitational instability. In- 
finite series of terms in standard Eulerian perturbation theory are resummed as a result of our starting from a 
Lagrangian description of perturbations. Delicate nonlinear effects on scales of the baryon acoustic oscillations 
are more accurately described by our method than the standard one. Our approach differs from other resum- 
mation techniques recently proposed, such as the renormalized perturbation theory, etc., in that we use simple 
techniques and thus resulting equations are undemanding to evaluate, and in that our approach is capable of 
quantifying the nonlinear effects in redshift space. The power spectrum and correlation function of our ap- 
proach are in good agreement with numerical simulations in literature on scales of baryon acoustic oscillations. 
Especially, nonlinear effects on the baryon acoustic peak of the correlation function are accurately described 
both in real space and in redshift space. Our approach provides a unique opportunity to analytically investigate 
the nonlinear effects on baryon acoustic scales in observable redshift space, which is requisite in constraining 
the nature of dark energy, the curvature of the Universe, etc., by redshift surveys. 
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I. INTRODUCTION 

Recent progress in cosmology is largely stimulated by pre- 
cise measurements of the Universe. The redshift survey of 
galaxies is one of the most important methods in distinguish- 
ing cosmological models. The baryon acoustic oscillations 
(BAOs) in the early universe Ol], which were first detected in 
the cosmic microwave background (CMB) anisotropy power 
spectrum |2], were also detected in the large-scale structure 
probed by redshift surveys IH. The expansion history and 
curvature of the Universe can be geometrically measured by 
galaxy clustering in redshift space Slst], and the BAO signa- 
ture in the large-scale structure yields a standard ruler 1^ for 
this purpose. Therefore, accurately measuring the BAO scales 
in a statistical quantity such as the power spectrum or corre- 
lation function of galaxies, one can constrain the structure of 
the Universe, such as the nature of dark energy [%. 

To utilize the BAOs as a standard ruler in the Universe, a 
theoretically precise description of the BAOs is of the utmost 
importance. The physics of the BAOs at sufficiently large red- 
shift is well described by linear theory. Only the amplitude of 
the spatial pattern of clustering is increased, and the character- 
istic scale of the BAOs is unaltered in this regime. However, 
the nonlinear dynamics play an important role in a redshift 
range where we can measure the large-scale structure. Even 
though the BAO scales ^ 100/j 'Mpc are large, the nonlin- 
earity deforms the BAOpattern in the power spectrum and 
correlation function lH |1 [T^, [H [H [Tj] . Nonlinear defor- 
mations of the power spectrum are described by the standard 
perturbation theory (SPT) |ll[ll[ll[ll[Il[ll|2^ of higher 
orders as long as the density fluctuations on BAO scales are 
still in the linear regime at sufficiently large redshift ll2ll I22I1 . 
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The last condition is not necessarily satisfied by realistic sur- 
veys. Moreover, when the cold dark matter (CDM) is present, 
the higher-order SPT does not give any sensible prediction for 
the correlation function because unphysical behavior in small 
wavelength limit prevents Fourier transform from converging. 

Recently, several other ap proa ches beyond the framework 
of SPT have been proposed Sllllllllilllllillim. 
Among others, Scoccimarro's reformulation of fluid equations 
using the propagator, the vertex, and a source ll23ll provides a 
way to use standard tools of field theory, yielding the renor- 
malized perturbation theory (RPT) [ZSl \2G, 27], the large- 

expansions ll28ll . the renormalization group approach ll29tl . 
and the closure theory |'30] in the context of gravitational in- 
stability theory. In these newly developed approaches, the 
standard perturbative expansion is reorganized and partially 
resummed in various ways. It is still difficult to exactly cal- 
culate the power spectrum by these approaches, and different 
levels of approximations and ansatz are employed. The pre- 
dictions of the power spectrum calculated by the 2-loop RPT 
and renormalization group method are shown to agree with 
measurements of the A^-body simulations on BAO scales even 
at low redshift where the 1-loop SPT breaks down [ 27, 290. 
It is still nontrivial why these approaches work well despite 
their approximate treatments and use of ansatz. 

So far the predictions of the above nonlinear theories, such 
as SPT, RPT, and their variants, have been confined in real 
space. Since the large-scale structure is actually measured 
in redshift space, it is mandatory to make predictions in red- 
shift space for applying the nonlinear theories to realworld 
observations. The redshift-space distortion effects in the lin- 
ear regime are anaMically given by first-order SPT (i.e., lin- 
ear theory) SIM Hi- However, investigations of nonlinear 
redshift-space distortions by higher-order SPT are quite lim- 
ited in literature Issi l34ll , apart from numerical simulations 
and nonlinear modelings ll 111 [H [H IH HI . 

In this paper, we consider resummations of the perturbative 
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expansion by a different approach from RPT and its variants. 
Instead of relying on Scoccimarro's reformulation, we start 
from a Lagrangian picture of density perturbations. Why do 
we have to add yet another approach even though there are 
akeady enough? There are a number of reasons. First, our 
approach is simpler than other resummation methods. While 
it is straightforward to calculate the power spectrum in our ap- 
proach, other resummation techniques are so complicated that 
approximate treatments and/or ansatz are necessary in deriv- 
ing the final power spectrum. Second, as a result of that, our 
analytic expression of the resulting power spectrum is no more 
complicated than the SPT expression, and numerical evalua- 
tions are fast performed. Third, and most importantly, our 
approach is capable of evaluating the nonlinear power spec- 
trum in redshift space, which is beyond the scope of other 
resummation methods in literature. Fourth, unlike in SPT, the 
power spectrum in our approach can be Fourier transformed 
to predict the nonlinear effects in the correlation function. 

An important ingredient of our approach is the perturbation 
theory in Lagrangian space, which is known as the Lagrangian 
perturbation theory (LPT) fi^. The first-order LPT corre- 
sponds to the classic Zel'dovich approximation jj/ll. Since 
the observable power spectrum is given in Eulerian space, pre- 
dicting the power spectrum by LPT is not a trivial task. We 
define a method of calculating the power spectrum from the 
LPT, partially expanding the Lagrangian variables in Eulerian 
space. The resulting expression contains an infinite series of 
SPT terms; that is, the SPT terms are automatically resummed 
in our formulation. The underlying motivation of this for- 
mulation is shared with RPT. Readers are recommended to 
go through an excellent introduction in Ref. ll25ll using the 
Zel'dovich approximation, which equally applies to our for- 
mulation since our formalism and RPT become equivalent in 
a case of Zel'dovich approximation. 

The rest of this paper is organized as follows. In Sec. [Ill 
our basic formalism is described. A brief review of LPT is 
also included. In Sec. |III1 analytic results of the power spec- 
trum and correlation function in real space are presented and 
compared with simulation results in the literature. The corre- 
sponding results in redshift space are derived in Sec. |IV] The 
conclusions are given in Sec.lVl 



n. RESUMMING PERTURBATIONS VIA THE 
LAGRANGIAN PICTURE 

A. The power spectrum from the Lagrangian variable 



Assuming the initial density field is sufficiently uniform, 
the Eulerian density field p{x) satisfies the continuity relation, 
p(x)d^x - pd^q, where p is the mean density in comoving 
coordinates. This relation is equivalent to the following equa- 
tion: 
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where 6 - plp-l'xs, the density contrast. The power spectrum 
P(k) is defined by a relation 
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5(fc) 



J d^xe-""-- 



■■6{x) 



(3) 



(4) 



is the Fourier transform of the density contrast. The power 
spectrum in real space depends only on the magnitude of the 
wave vector, k - \k\. In redshift space, however, a directional 
dependence emerges. Using Eqs. ©-(Hli, the power spectrum 
can be expressed by a displacement field as 13811 
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where q - q\ - qi- Because of the homogeneity, the en- 
semble average in the integrand of the right-hand side (RHS) 
only depends on the separation q. The expression of Eq. (|5]l is 
a general relation between the power spectrum and displace- 
ment fields. 

We use the cumulant expansion theorem lf39ll 



{e-'^} = exp 



i-iT 
Nl 
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in Eq. (|5]l, where {X^)c denotes a cumulant of a random vari- 
able X 1 20]. In the resulting equation, we have a following 
quantity: 



{{k ■ - ^(q2)]f)^ 

= [i+(-if]([fc-^(0)f)^ 



N-l / > 



(7) 



In a Lagrangian picture, a fundamental variable to repre- 
sent perturbations is a displacement field 'ifiq, t), which maps 
a fluid element from initial Lagrangian coordinates q to the 
Eulerian coordinates a; at a time f: 



x{q, t)^q + t). 



(1) 



We will drop the argument t in the following for notational 
simplicity, and readers should keep in mind that the dynamical 
variables always depend on time. 



where the binomial theorem is applied. In this equation, only 
N >2 cases are survived, because {^} = for parity symme- 
try. Accordingly, we obtain a general expression. 
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where summation over repeated indices is assumed as usual, 
and 



X (qi) ■ ■ ■ (92) ■ ■ ■ *l',„)(g2))^ 



(9) 
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are defined. In the RHS of Eq. ( fTOl l. all the indices are 
symmetrized over 

The quantity A^^* of Eq. ^ is given by a cumulant of a 
displacement vector at a single position. The quantity B*'^' of 
Eq. ( fTOb is given by cumulants among two displacement vec- 
tors which are separated by a Lagrangian distance |q|. Corre- 
lations among components of separated displacement vectors 
are expected to be small when the separations are large. When 
we are interested in clustering on large scales, the magnitude 
of wave vector k is small, and large separations of |<7| in Eq. (|5]l 
efficiently contribute to the integral. Therefore, we expect the 
cumulants of B'^'^' are smaller than the cumulants of A"^' on 
large scales. 

With the above consideration, we propose to expand only 
the B^^^ terms in Eq. dS]! as small quantities and keep the A*'^' 
terms in the exponent. If the A^'^^ terms are expanded as well, 
this equation turns out to give equivalent results of SPT, which 
is formulated in Eulerian space, as shown below. Therefore, 
the exponential prefactor contains infinitely higher-order per- 
turbations in terms of SPT. Thus our method provides a way 
to resum the infinite series of perturbations in the SPT. This 
point is our first key technique of this paper. The second tech- 
nique is to evaluate the cumulants of the displacement field by 
the LPT, which is addressed in the following. Before that, it 
is convenient to represent the cumulants of the displacement 
field in Fourier space. 



B. Poly spectra of the displacement field 

To evaluate Eq. (jS), it is useful to define the polyspectra 
C,,...,„ of the displacement field by 

(t„(pi)--->I',„(p^))^ 

^{2nf6\pi+---+pN){-if-^Ci,..,,{p,,...,pN), (11) 

where 
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is the Fourier transform of the displacement field. The re- 
lation pi + ■ ■ ■ + pm - is always satisfied because of the 
translational invariance. For N - 2, we employ a notation, 
Cijip) - Cijip, -p), for simplicity. The factors (-i)^^^ in the 
RHS of Eq. ( fTTT i are there to ensure that the polyspectra C,, ... 
are real numbers. In fact, one easily finds the complex conju- 
gate satisfies 



(^^.(Pl) ■ ■ ■ *,,(Piv))* = (-1)^ (Pi) ■ ■ ■ "^iAPN))^ , 



which is an equivalent condition that the polyspectra are real 
numbers, C*...,^^ = Ci^...i^. It is also important to note the 
polyspectra of the displacement field satisfy a parity relation, 

C,-,...,-„(-pi, . . . , -Pn) = (-l)'^C,-,...,-„(p,, . . .,Pn), (14) 

which is explicitly shown by Eqs. ( fTTT i and ( fT2b . 

Using the polyspectra defined above, Eqs. ^ and ( fTOl ) re- 
duce to 
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C. The Lagrangian perturbation theory 

For a Newtonian pressureless self-gravitating fluid embed- 
ded in an expanding universe, the displacement field is gov- 
erned by the equation of motion. 



d^^ c/* 
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where H - a/ a is the time-dependent Hubble parameter, and 
V I is a derivative with respect to the Eulerian coordinate, x = 
q + '4'(q). The gravitational potential (p is determined by the 
Poisson equation. 



Vicpix) = 4nGpa'^6{x) 



(18) 



where the density contrast is related to the displacement field 
by Eq. (O or, equivalently. 



6{x)^{det[6u + 'i'ij{q)]Y' - I. 
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We use the LPT f36'] to evaluate the polyspectra of the dis- 
placement field. In LPT, the displacement field is expanded 
by a perturbative series: 



(20) 



where ^^-^^ has the order of ('4'^'')'*', which is considered 
a small quantity (properly speaking, dimensionless quantity 
di^j is considered to be small). It is quite common in the LPT 
that the velocity field is assumed to be irrotational. 



V. X I) = 0. 



(21) 
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This condition is compatible with the dynamical Eq. ( [TtI i and 
restricts one to considering a subclass of general solutions. 
The main motivation of this restriction is that the vortical 
motions linearly decay in Eulerian perturbation theory It2dll 
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and thus are negligible in the quasilinear regime. We assume 
this condition throughout this paper The longitudinal part of 
Eq. (TJl is an independent equation to be solved perturbatively 
in LPT; the same is true for Eqs. ( fTSl l. ( fT9] l. and (1211 1. Each 
perturbative term of the displacement field in Fourier space is 
generally represented as 



iD" r cPp 
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■6o(p„), (22) 



where (Jo - 5^''(f = fo) is the linear density contrast at the 
present time to, D{t) is a linear growth rate normalized by 
D{t()) = 1. The imaginary unit is there to ensure that the per- 
turbative kernels X*"* are real vectors. This is straightforward 
to show by considering the parity transformation of both sides 
of Eq. (I22I 1. The perturbative kernels do not depend on time 
in an Einstein-de Sitter model. In general cosmology, they 
weakly depend on time and cosmological parameters |40]. It 
is still a good approximation that the perturbative kernels in 
arbitrary cosmology are replaced by those in the Einstein- 
de Sitter model [20]. We adopt this approximation throughout 
this paper for simplicity, although including such dependence 
is straightforward. The perturbative kernels in LPT up to third 
order are given by ll4lll 
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where k - pi + ■ ■ ■ + p„ for L^"\ and a vector T represents a 
transverse part whose expression is not needed in the follow- 
ing application. It is useful to symmetrize the kernel i^^^' in 
terms of their arguments: 

i^^'(Pi,P2,P3) = ^[i'^"^(pi,P2,P3) + perm.]. (26) 



D. One-loop approximation to tlie resummed power spectrum 

The polyspectra of the displacement field in LPT can be 
evaluated in almost the same way as in the case of density 
fields in Eulerian perturbation theory. Substituting Eqs. ( fT2b . 
( I20I 1. and ( I22I 1 into Eq. ( fTTT i. the perturbative expansion of the 
polyspectra is calculated. We assume the initial density field 



is a random Gaussian field, in which case the cumulants of 
the linear density field are completely expressible by a linear 
power spectrum, Ph(k) = D^Po(k). Up to second order in 
Ph{k), we obtain 



Cijip) = c[)'>(p) + c)f'{p) + c:/>(p) + cy/'ip) 
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where we define quantities c|j"'\ C^™'' by 

{^'f\p)¥p\p')l = (2n)'6\p + p')d;;'"\p), 

{¥;'\piWf\P2)'^f(P3)l 

= (27Tf6\pi+P2+P3)Cf.f(puP2,P3), 

which are exphcitly given by using LPT kernels as 
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and 



d:,"'(pi,P2,P3) = C*^"'(p3,Pl,P2) 
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= -L)"{pi)Lf{p2)q'{puP2)Papy)PaP2). (34) 

The polyspectra of order four or higher do not contribute up 
to second order in P\^{k). In random Gaussian fields, it can 
generally be shown that A^-polyspectrum has the order - 1 
in P^ik), for the connectedness of the cumulants. In SPT, dia- 
grammatic techniques prove to be quite useful llT7ll42ll . Simi- 
larly, we find diagrammatic representations of the above quan- 
tities of Eqs. OD-Oll, which are given in Fig.[T] 

At this point, Eq. (O is expanded by LPT by using 
Eqs. ( fTSl ), ( fT6l ), and (I3n)-(l34l). As promised, only quantities 
Z?('^' are expanded from the exponent in Eq. ([8j and keep A^'^' 
in the exponent. The result is 

/j3 
^4''>(p) 

+ d'f '(fc, -p, p-k) + cf^ '\k, -p, p - k)] 

+ \kikjhk, j ^cf/\p)Cf,'\k-p)^. (35) 
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pansion of Eq. ( l35T l is therefore a loop expansion, and all the 

1- loop diagrams are exhausted in this equation. In the expo- 
nent of the prefactor, we could consider the higher-order terms 
such as integrals of C^^^\ Cfj^K These terms correspond to the 

2- loop diagrams of Fig. [3] Although those terms are second 
order in P^ik), their contributions to the total power spectrum 
are third order since they are multiplied by first- or higher- 
order terms in Eq. (l35b . Therefore, we do not include these 
2-loop contributions and truncate all the multiloop contribu- 
tions. 



FIG. 1: Diagrammatic representations of C.j (p) (a), C.j (p) (b), 
C^^ip) (c), and C<;'''(p,,p2,p3) (d). 



Except for the linear-theory contribution, cfj^\k), all the 
polyspectra in the above terms are accompanied by integra- 
tions over wave vectors. Accordingly, these terms correspond 
to the 1-loop diagrams in Fig. |2l The bubble diagram (a) of 

Fig.|2]corresponds to the exponent of the first factor. The tree 

-,(11) 



diagram (b) corresponds to the linear contribution C.j , and 
loop diagrams (c)-(e) correspond to C|^^', C|'^', and C;; re- 

'J U 'J 

spectively. The loop diagrams (f)-(h) correspond to integrals 
of Cfjl^\ d]f*, and Cfjl^\ respectively. The 
gram i) corresponds to the last integral in Eq. 



III. CLUSTERING IN REAL SPACE 

A. The power spectrum 

It is straightforward to calculate Eq. ( |35] | using the expres- 
sion of the polyspectra of the displacement field, Eqs. OTTl- 
( |34| ). and perturbation kernels, Eqs. (l23Tl-(l26]l. Contractions 
with spatial indices in Eq. ( |35] | can be taken before evaluat- 
ing the momentum integrations. The transverse part in the 
third-order kernel, Eq. dZSl ), vanishes in course of calculation, 
because of the inner product with vector k. As a result, the 
calculation is quite similar to that in SPT, where some of the 
angular integrations can be analytically done 1 18]. The result 
is 



P(k) = exp 

X [pdk) + 
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Comparing this result with that of 1-loop SPT flS*!, one im- 
mediately notices the similarity of the equations. The only 
differences are the existence of the exponential prefactor and 
the number +10 of the second term in the last square bracket, 
which is -158 in SPT instead. Thus, Eq. (|36] | is equivalently 
represented as 



P(k) = exp 



P^k) + Psp7(^) + 



where Pgp!J°''(fc) is a pure 1-loop contribution (without a linear 
contribution) to the power spectrum in SPT. It is obvious from 
this expression that expanding the exponential prefactor ex- 
actly reproduces the SPT result at the 1-loop level. The expo- 



J dpPdp) 

k^ r 

^PUk) J dpPdp) 
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nential prefactor contains information from an infinite series 
of higher-order perturbations in terms of SPT. In terms of dia- 
grams of SPT 111711 . the exponential prefactor comes from the 
diagrams in which two points are connected by a single inter- 
nal line with P^ik) but vertices are dressed by loops. However, 
those diagrams are not completely resummed, and parts of the 
resulting terms are resummed. Therefore, the resummation of 
this paper corresponds to a partial vertex renormalization in 
terms of the SPT. 

On the other hand, the resulting expression has a certain 
similarity with RPT |25, 26], because of the exponential pref- 
actor In fact, a square root of this factor coincides with the 
nonlinear propagator in a large-fc limit of RPT. In its formula- 
tion, the nonlinear power spectrum is divided into two parts. 
The first part is linearly dependent on the initial power spec- 
trum with a coefficient given by the square of the propagator 



j p-k 



p — k 



0.01 



0.1 

k [hMpc-i] 



FIG. 4: Comparison of power spectra by different approximations 
at redshifts z = 0, 1,2,4 (from upper to lower lines). Black (solid) 
line: this work; red (dotted) line: linear theory; green (dashed) line: 
1-loop SPT. Two nonlinear scales ^nl, ^nl/S in each redshift are in- 
dicated by arrows. 



k - p \ / k — p 

- - * * 
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i) 



FIG. 2: Diagrams for our resummed power spectrum of Eq. l |35t . 




b) 



FIG. 3: Two-loop bubble diagrams (see text). 



The second part corresponds to contributions from mode cou- 
plings among different scales. RPT uses loop expansions to 
evaluate the mode-coupling contributions. Our expression has 
a similar structure, although the methods are quite different. 
One possible reason for this similarity is that the RPT propa- 
gator in the Zel'dovich approximation exactly coincides with 
the square root of the exponential prefactor of Eq. ( [36l ) above. 
In the Zel'dovich approximation, only c!j'' is present for cu- 
mulants of the displacement field. As a result, loop corrections 
come only from diagrams of type i) in Fig.|2]and its multiloop 



extensions («-loop diagrams with «-t- 1 internal lines connected 
to the two external vertices), in which the vertex renormaliza- 
tion cannot appear. As a result, RPT and our approach turn 
out to be equivalent in the Zel'dovich approximation (see the 
introduction of Ref. i25ll ). 

Having an exponential damping factor, the power spectrum 
in our approximations does not have power at large wave 
numbers. This behavior is not physical in the fully nonlin- 
ear regime, since the nonlinear evolution enhances the power 
spectrum on small scales. Our approximations only apply to 
the quasilinear regime on large scales. An advantage of this 
property is that the integration of the power spectrum is stable 
so that other statistical quantities may be obtained, such as the 
correlation function, on quasilinear scales. 

On the other hand, our approximations for the power spec- 
trum itself fail when the exponential damping factor is effi- 
cient. It is useful to have an estimate for the validity range 
of wave numbers in our approximations. Since the damping 
scale is given by 



dpPdp) 



-1/2 



(38) 



a possible criterion for the validity range is A: < ak^^i^, where 
a should be determined according to the required accuracy. 
We will see below that a = 0.5 corresponds to the accuracy 
within a few percent, compared to the A^-body simulations. 
Throughout the rest of this paper, we consider k < A:nl/2 as 
the validity range in our approximations. 

In Fig.m the power spectra of our result, and those of linear 
theory and 1-Ioop SPT, are compared. The linear power spec- 
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trum Pdk) is calculated from the output of CMBFAST 
at z = 0. We adopt a cosmological model with parameters 
fin, = 0.27, Qa = 0.73, Qb = 0.046, /z = 0.72, = 1 and 
erg = 0.9. The damping scale A:nl and the claimed validity 
scale A;nl/2 are indicated in the figure. It is clearly seen that 
our power spectra are expnentially damped on small scales 
k > A;nl- In this log-log plot, differences among the three 
approximations in our claimed validity range k < ^nl/2 are 
hardly seen. Our approach is advantageous when a very accu- 
rate estimation on large scales is crucial. The determination 
of the BAO scale really requires such accuracy. 

In Fig. |5] detailed features in nonlinear power spectra are 
shown on BAO scales. To enhance the baryon oscillations, 
each power spectrum is divided by a smoothed, no-wiggle 
linear power spectrum Pnw(^) |44]. The same cosmological 
parameters are assumed as in Fig.|4] The validity ranges are 
indicated in the figure. There are differences among three ap- 
proximations in our validity ranges. 

In Fig. |6l our result of the power spectrum in real space 
is compared with the A^-body simulations of Ref. JH] using 
the same parameters as Seo and Eisenstein's. In this case, 
ris - 0.99 and other cosmological parameters are the same 
as above. The linear power spectrum is calculated from the 
output of CMBFAST at z = 49 to match the input power spec- 
trum of the simulations. Since the no-wiggle power spectrum 
Pnviik) matches to the power spectrum at z = 0, the baseline is 
slightly larger than unity, because of nongravitational effects 
in the transfer function. To reduce the effects of cosmic vari- 
ance in A^-body simulations, we calculate the differences (in 
logarithmic scales) between the output power spectra and lin- 
early evolved initial power spectra at z = 49, both measured in 
simulations. The resulting differences are added to the linear 
theory prediction in the figure. In this way, initial deviations 
of the power spectrum from the theoretical power spectrum 
are subtracted (similar presentations are seen in Ref. |9]). The 
actual output power spectrum of simulations is noisier than 
that in the figure especially on larger scales. As described 
in Ref. fs*], the Savitzky-Golay filtering is applied when the 
power spectrum in each A^-body sample is evaluated. The shot 
noise is not corrected for the A^-body samples, since the initial 
positions of particles are on regular grids, and the naive shot 
noise Psn(^) - i/n, where n is the number density of parti- 
cles, does not apply at a large redshift in the simulation Il45[| . 
It is not obvious how to correct for the shot noise in such a sit- 
uation. The samples with z - 0.3 and z = 1 are hardly affected 
by the shot noise because the signal amplitudes are large, but 
contributions from the naive shot noise in the z = 3 sample are 
over 1% for k > 0.094/1 Mpc"'. Therefore, the z = 3 sample 
could be more or less affected by nontrivial shot noise in that 
range, although that is smaller than the naive shot noise. 

In all cases of the three redshifts z = 0.3, 1, 3, the ranges 
where our results and simulations agree are broader than those 
of 1-loop SPT. Overall, the 1-loop SPT overestimates the 
power spectrum on relatively small scales, while our result 
is always better on the same scales. For example, the height 
of the peak at k ^ 0.07 h Mpc ' of the z = 0.3 sample is ex- 
plained by our work, while the 1-loop SPT overestimates the 
height. In a range k ^ 0.09-0.13 /i Mpc ' of the z = 1 sam- 



ple, A^-body data are (accidentally) closer to the linear theory 
rather than the 1-loop SPT. This tendency is also explained by 
our work. In the z = 3 sample, our result matches the A^-body 
result in the range k < 0.2 h Mpc while linear theory and 1- 
loop SPT do not (under the condition that the nontrivial shot 
noise is negligible as noticed above). The simulations and 
our results agree within a few percent in our claimed validity 
range, k < /:nl/2. 

B. The correlation function 

The correlation function ^(r) is given by a Fourier transform 
of the power spectrum. After angular integrations, they are 
related by 

X°° k^dk 
—j,(kr)P(k) (39) 

where joix) = jc ' sinx is the spherical Bessel function of 
zeroth order. 

In the 1-loop SPT, the power spectrum badly behaves on 
small scales to perform this integral. The 1-loop contribu- 
tion to the power spectrum in SPT consists of two terms, 
Pjp°°''(fe) - Pnik) + Pn(k). Each term asymptotically be- 
haves as III El 

k^p,(k) r°° 

Piiik) ^ -Pnik) ^ dpPi^(p\ {k^^), (40) 

and the net power spectrum is much smaller in magnitude than 
each term. On mildly small scales, the residuals also behave 
as oc k^Pi^(k) for a CDM-like power spectrum. Therefore, the 
power spectrum on small scales contributes to the integral in 
Eq. ( [39] l even if the separation r is large. In a small-scale limit, 
k oo, the leading contributions of the two terms are exactly 
canceled out, and the net contribution behaves as lITsilTsIl 

p;j!°°P(A:) (X PL(fc) dpp^Pdp), (k^<x,). (41) 
Jo 

For a CDM-like power spectrum, d In Pi^(k)/ dink > -3, and 
the integral in Eq. ( |4T] ) diverges. As a result, one cannot eval- 
uate the correlation function by the 1-loop SPT, because the 
integral of Eq. ( [39] l does not converge. This divergence obvi- 
ously reflects the fact that the 1-loop SPT is invalid on small 
scales. Since the correlation function is an observable quan- 
tity, this is an undesirable property in SPT. Such a problem 
does not occur in our expression of Eq. ( [36l ): the power spec- 
trum is exponentially suppressed in the small-scale limit. This 
behavior is not physical on sufficiently small scales, since our 
expression is applicable only to the weakly nonlinear regime. 
However, small-scale behavior in the power spectrum should 
not be relevant to the correlation function on large scales, 
since dominant contributions of the power spectrum to the cor- 
relation function of Eq. ( [39l ) come from scales ~ r ' for a 
CDM-like power spectrum. 

Since the exponential prefactor in Eq. ( [36] l acts as a Gaus- 
sian damping function, the resulting correlation function has a 
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FIG. 5: Nonlinear evolution of the baryon acoustic oscillations in real space for various redshifts, z = (top left), 0.5 (center left), 1 (bottom 
left), 2 (top right), 3 (center right), 4 (bottom right). Each power spectrum is divided by a smoothed, no-wiggle linear power spectrum P„„(k) 
I44I1 . Black (solid) line: this work; red (dotted) line: linear theory; green (dashed) line: 1-loop SPT. The validity range k < A;nl/2, where our 
result is expected to be accurate within a few percent, is indicated by arrows. 



form of convolution by a Gaussian smoothing function. This 
smoothing effect by nonlinearity is already known, which 
smears the BAO peak in the correlation function 01311 . The 
physical origin of this nonlinear smearing is the random mo- 
tion of mass elements, because the dominant contribution of 
this smearing factor comes from (l'4'P) as obviously seen by 
Eqs. ([8]l and (|9j. It is worth noting that this nonlinear smearing 
effect is not captured by linear theory and higher-order SPT. 

In Fig. |7] the nonlinear correlation function is plotted on 
BAO scales. The same cosmological parameters as in Fig. |5] 
are adopted. Nonlinear smearing effects on the baryon acous- 
tic peak are described by our theory. In Fig. |7] Gaussian- 
filtered linear theory is overplotted. The filtering length is 
given by the one that the exponential prefactor in Eq. ( [36] l 
implies: 

nl/2 



dpPUp) 



(42) 



This simple model approximately explains the nonlinear 



smearing effects. In detail, this simple model overestimates 
the small-scale clustering and height of the trough, and under- 
estimates the peak height. 

In Fig. |8] our result of the coiTelation function in real 
space at z - 0.3 is compared with the A^-body simulations 
of Refs. 1^ [4^ using the same parameters as used in these 
references. As in the case of the power spectrum, we compute 
differences (in linear scales) between the output correlation 
function and the linearly evolved initial correlation function 
at z = 49 in the simulation, and the resulting differences are 
added to the linear-theory prediction in the figure. If we di- 
rectly plot the measured coiTelation function, the amplitudes 
on r ^ 70-100 /j^'Mpc and 1 15-135 /i^'Mpc are slightly 
(about 10"^) larger This tendency is already seen even in the 
initial correlation function at z = 49, which can be identified 
to the finiteness effect of the simulations. Therefore, we only 
measure the nonlinear deviations from the linear theory in the 
simulation. 

The agreement between the simulation and our result is 
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FIG. 6: Comparison of the power spectra to the W-body simulations 
of Ref. 1 8] in real space at redshifts z = 0.3, 1 , and 3 from top to bot- 
tom, respectively. Open circles: A'-body simulations; Black (solid) 
line: this work; red (dotted) line: linear theory; green (dashed) line: 
1-loop SPT. Only nonlinear deviations from the linear growth is mea- 
sured in A'-body simulations to reduce cosmic variances. Estimates 
of validity range, k < A:nl/2 are indicated by arrows. In the z = 3 
sample, the range where contributions from the naive shot noise ex- 
ceed 1% is also indicated, althogh the true shot noise could be smaller 
because of the regularity of particle positions at such a large redshift 
in the simulation. 



quite good. There is not any fitting parameter at all. In the 
corresponding power spectrum, the z - 0.3 plot in Fig. |6l 
the difference between our result and that of 1-loop SPT may 
seem to be small. In the correlation function, however, the dif- 
ference is really big, since the 1-loop SPT cannot predict the 
correlation function, while our result explains the nonlinear 
effects on BAO scales very precisely. 



IV. CLUSTERING IN REDSHIFT SPACE 
A. Mapping from real space to redshift space 

Our approach above can be extensible to the calculation in 
redshift space. In redshift surveys, the position of each galaxy 
is measured by the redshift. In a completely homogeneous 
universe, the comoving distance-redshift relation is given by 
a standard formula for the Friedmann-Lemaitre universe: 



x(z) 



cz 



(43) 



where H{z) is the time-dependent Hubble parameter. In re- 
ality, however, the Universe is locally inhomogeneous, and 
the comoving distance to a galaxy of redshift z is not exactly 
given by Eq. ( |43] ). For the galaxy redshifts, the leading source 
of the deviation is the Doppler shift caused by the peculiar 
velocities of the galaxy. The observed redshift Zobs with the 
effect of the peculiar velocity can be calculated by solving the 
geodesic equation. The comoving distance in redshift space 
is defined by s - Ji:(Zobs), where, in the RHS, the redshift- 
distance relation of Eq. ( l43T l is adopted even though we live 
in an inhomogeneous universe. When the peculiar velocities 
are nonrelativistic, the comoving distance in real space, x, and 
that in redshift space s is related hy s - x + v^JaH |47], where 
V, - V ■ z \s ?L line-of-sight component of the peculiar velocity 
V. Throughout this paper, we work in "distant-observer" or 
"plane-parallel" approximation, where the line-of-sight direc- 
tion is fixed in space and denoted by z. We ignore the contri- 
bution from the peculiar velocity of the observer, which does 
not contribute to the following calculation of the power spec- 
trum. When the curvature scale of the Universe is much larger 
than clustering scales we are interested in, we can safely use 
local Cartesian coordinates in the volume of galaxy surveys. 
Thus the relation between the position in real space x and that 
in redshift space s is given by 



z ■ V 

S - X -i Z, 

aH 



(44) 



where v - ax. Therefore, using Eq. ([T), the displacement 
field in redshift space is given by 



z ■ # 
H 



(45) 



A similar equation has been applied to the analysis of the 
Zel'dovich approximation |48]. 

In the time-independent approximation to the perturbative 
kernels, 'i'*"' oc D", and therefore, the time derivative of the 
displacement field is simply 



(46) 



where / = dlnD/dlna - {HD)^^D is the logarithmic deriva- 
tive of the linear growth rate. Thus, the displacement field of 
each perturbation order in redshift space is given by 



(47) 
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FIG. 7: Nonlinear evolution of the baryon acoustic peak in real-space correlation function for various redshifts, z = (top left), 0.5 (top right), 
1 (bottom left), 2 (bottom right). Black (solid) line: this work; red (dotted) line: linear theory; green (dashed) line: Gaussian-filtered linear 
theory. The 1-loop SPT cannot predict the correlation function. 
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FIG. 8: Comparison of the correlation functions to the A^-body sim- 
ulations of Refs. 1 8, 46] in real space. Open circles: A'-body results; 
Black (solid) line: this work; red (dotted) line: linear theory. Only 
nonlinear deviations from the linear growth are measured in A'-body 
simulations to reduce finiteness effects. 



with which Eq. (|47| reduces to ^'f"^ = R^f¥"\ or in a vector 

' U J 

notation, = Therefore, the perturbative kernels 

in redshift space are given simply by 



' s(n) 



(«) r (") 



(49) 



Thus, our calculation in real space can be simply generalized 
to that in redshift space by using the redshift-space perturba- 
tive kernels in Eqs. (I3ni-(l34b. while the form of Eq. dSST l is 
unchanged in redshift space, provided that the cumulants of 
the displacement field are evaluated in redshift space. The 
mappings of the order-by-order cumulants are quite simple: 



C 



s(«m) 



- ^ik ^jl ^kl 



^s(nin2K3) _ j^lm) j^(n2) 



'Ijl l2]2 hj3 



iihh ' 



(50) 
(51) 



B. The power spectrum in redshift space 



which is just a linear mapping of the displacement vector of 
each order This linear transformation is characterized by a 

iii) 

redshift-space distortion tensor R.. for nth order perturba- 
tions, defined by 



The calculation of Eq. ( [35] ) in redshift space is more com- 
plicated than in real space, since the anisotropy is introduced 
by the redshift-space distortion tensors. Accordingly, the mo- 
mentum integrations should be evaluated before taking in- 
ner products with the vector fc. In redshift space, however, 
the cumulants C*"'"* and in Eq. (O should be re- 



(48) placed by redshift-space counterparts, Cjj 



and C^<"'"^"-'' of 
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Eqs. ( [soi l and ( fSTT ). Since the redshift distortion tensors R^"^ 
are anisotropic, the integrands do not reduce to scalar func- 
tions. Therefore we need to explicitly evaluate integrals of 
tensors. The results are given in Appendix |A] Eqs. ( |A9t - 
(lAl 11 1 and (IA13b -( lA15l l. Using those equations, it is straight- 
forward to obtain an expression of the power spectrum P^ik) 
in redshift space. There are many terms in the result, which 
are arranged as 



Ps(fe) = exp{-fc2[l+/(/ + 2)/z2]A} 



(52) 



where = z ■ k/k is the direction cosine of the wave vector k 
with respect to the line of sight. 



dpPdp), 



(53) 



and 
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(55) 
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57^ 51 
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(57) 


^23 


= -y2i-322 + ^24- 


7 ' 




(58) 


£24 








(59) 


£33 


3 27 
= y2i + y22+223 - 


^24 + ^^1 + 


12 

—R2, 


(60) 


£34 


3 3 1 
= -321-522 + 424, 






(61) 


£44 


= ^21 + ^22 + ^23- 


^24, 




(62) 



and all the other E„,„ which are not listed above are zero. The 
functions Q„{k) and R„{k) are given by Eqs. ( lAlb - dASl l in Ap- 
pendixlAl 

As in the case of real space, the SPT result should be re- 
produced by expanding the exponential prefactor in Eq. ( |52] |. 
A formal expression for the 1-loop power spectrum in red- 
shift space iyy\ contains 3 -dimensional integrations. Further 
reduction of the expression is possible and necessary to com- 
pare with our result. In Appendix |B] we derive a reduced 
expression for the power spectrum of 1-loop SPT in redshift 
space, which itself is a new result. Finally, the following rela- 
tion, which is similar to Eq. (|37] |. is confirmed: 

p,{k)^tx^[-k\\+f{f + 2)^?]A] 
x[{\+f^?fp,^(k) + p';^'^^{k) 

+ (1 + fulfil + fif + 2)i?]k^Pak)A} , (63) 



where P^'^°°''(k) is the 1-loop contribution to the redshift- 
space power spectrum in SPT given in Appendix |B] There- 
fore, when the exponential prefactor is expanded, Eq. ( |52| | re- 
duces to the result of SPT. 

The linear term in Eq. (|52] | corresponds to the Kaiser for- 
mula ll3lll of the linear perturbation theory. The appearance of 
the exponential suppression prefactor is characteristic in our 
formalism. This factor reflects the nonlinear smearing effect 
of the velocity field along the lines of sight ll49ll . such as the 
fingers-of-God effect by large-scale random motions. The lin- 
ear theory and higher-order SPT do not have such a suppres- 
sion factor along the lines of sight. 

This form of the suppression factor exactly matches a phe- 
nomenological model of the redshift-space power spectrum on 
large scales [13], which is given by (after correcting obvious 
typos) 



P(k) - exp 



-Uklal^kjal) 



(l+ffi^fPdk), (64) 



where kj_ and k\\ are wave vector components across and along 
the line of sight, and cr± and cry are modeled by 

0-^ = soD, cr|| = soD(l + /). (65) 

The length parameter so characterizes the scale of radial dis- 
placements. When all the 1-loop contributions of our Eq. ( l52b 
are omitted, the model of Eq. ( l64b is "derived" if we identify 
5() = (2A)'^^/D - (2Ao)^^^, where Ao is a linear-theory value 
of A at z = 0. In Ref. ifTsll . the parameter sq is measured 
from the A^-body simulation of Eisenstein, Seo, and White, 
and a value sq - 9.40/i^'Mpc is obtained for 100/i^'Mpc 
separations in their cosmology (we convert their quoted value 
according to different conventions of the growth factor). Us- 
ing the same set of cosmological parameters as theirs, we find 
Ao = 44.0/i"^Mpc^. The corresponding prediction from our 
formalism is sq = 9.38 /i^'Mpc, which is very close to their 
measured value. 

There is another simple model of the nonlinear power spec- 
trum in redshift space |35], which is given by 

PJk) ^ exp{-fkVcrl) 

X [Pss(k) + 2fii^Psg(k) + fii'Pge(k)] , (66) 

where cr^ is the same as our A of Eq. (53[ . and Pssik), Pseik), 
Peeik) are (cross) power spectra of density and velocity fields, 
calculated from nonlinear 1-loop SPT. Even though the den- 
sity and velocity (cross) correlations are taken into account, 
this model is no longer a result of exact 1-loop perturbation 
theory (see AppendixlBli. This simple model explains A^-body 
data to some extent but predicts a slightly larger power spec- 
trum than A^-body data on small scales IBSll . This model has a 
similar structure to our result but lacks many terms. 

The spherical average of Eq. ( |52] i is calculated by using the 
following integral: 



(-1)" 





erf(x'/2y 







(67) 
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where -y(a, x) is the lower incomplete gamma function, and erf 
is the error function normalized by erf(+oo) - 1. 

In Fig. |9] the spherically averaged nonlinear power spec- 
trum in redshift space is plotted on BAO scales. The same 
cosmological parameters as in Fig.|4]are adopted. Each power 
spectrum is normahzed by a no-wiggle power spectrum lliill 
in real space. According to the Kaiser effect, the amplitude is 
larger in higher redshifts. In the same plot, the simple model 
of Eq. (l66T l is overplotted. Both linear theory and 1-loop SPT 
are larger than our result and the simple model, since the for- 
mer two do not have a suppression factor which corresponds 
to nonlinear smearing effects along the lines of sight in red- 
shift space. The validity ranges k < kf^hf^ are indicated in the 
figure, where ^nl is given by Eq. (|38] |. 

In Fig. (To] a comparison of the power spectrum in red- 
shift space between our results and the A^-body simulations 
of Ref. |8j] is presented. As in real space, our results are better 
than 1-loop SPT results. The 1-loop SPT overestimates the 
power spectrum much worse in redshift space. On smaller 
scales, the simple model of Eq. ( |66] l and the linear theory 
seems to be good descriptions at a particular redshift. How- 
ever, that is just a coincidence because they do not describe 
the simulations at different redshifts. The shot noise is not 
corrected for the A^-body simulation, and the A^-body sample 
of z = 3 might be affected by the shot noise on small scales. 
The range where the naive shot noise is over 1 % is indicated 
in the figure. Since the particle positions in redshift space are 
more randomly distributed than in real space, the true shot 
noise might be larger than that in real space. Overall, our re- 
sults describe the A^-body results well at all redshifts. 



C. The correlation function in redshift space 

The correlation function ^^if) in redshift space is given by 

/tJ^k 
-—e'^^-p^ik). (68) 

When we consider the spherically averaged correlation func- 
tion, the relation of Eq. ( [39] l is applicable with the spherically 
averaged power spectrum. 

In Fig. [TT] the spherically averaged nonlinear correlation 
function in redshift space is plotted on BAO scales. The same 
cosmological parameters as in Fig. |5] are adopted. Since the 
effect of the suppression factor is larger in redshift space, 
the nonlinear smearing effects on the baryon peak are larger 
than in real space (cf. Fig. |7]). This is because the nonlin- 
ear redshift-space distortions additionally smear the clustering 
along the line of sight. As in Fig. |7] Gaussian-filtered linear 
theory is overplotted, which approximately explains the non- 
linear smearing effects. 

In Fig. [12] our result of the correlation function in redshift 
space at z = 0.3 is compared with the A^-body simulations 
of Refs. Isl, 143] using the same parameter as used in these 
references. The agreement is quite good also in redshift space, 
and our result explains the nonlinear effects on BAO scales 
very precisely even in redshift space. 



V. CONCLUSIONS 

In this paper, we have formulated a new resummation 
method of the cosmological perturbation theory, starting from 
the Lagrangian description of perturbations. The infinite se- 
ries of perturbations in SPT is naturally contained in the ex- 
ponential prefactor in our formalism, and expanding the expo- 
nential factor reproduces the SPT. We have presented 1-loop 
results of our formalism in this paper. Our method is particu- 
larly useful for studying nonlinear effects on BAO scales and 
always provides a better description of the A^-body results than 
SPT on large scales. 

Nonlinear corrections to the correlation function can be cal- 
culated by our formalism. This is a great advantage over SPT, 
since higher-order SPT cannot predict the correlation func- 
tion. Our 1-loop results are comparable to the 1-loop results 
of the RPT 1 27]. In the RPT, it is shown that taking into ac- 
count the 2-loop effects enlarges the applicable range toward 
small scales in the power spectrum. Since our formalism and 
RPT have common features in the resulting equations, we ex- 
pect that 2-loop extensions of our calculation similarly enlarge 
the applicable range of scales in the power spectrum. The de- 
scription of the BAO peak in the correlation function seems 
already sufficient by our 1-loop calculation. 

It is an important and original outcome that our formal- 
ism gives an analytical description of nonlinear effects in red- 
shift space. The large-scale structure is measured by redshift 
surveys of astronomical objects, such as galaxies, quasars, 
Lyman-alpha absorbers, 21 cm emission, etc. However, non- 
linear corrections to the power spectrum and correlation func- 
tion in redshift space by perturbation theory have not been suf- 
ficiently studied so far Our formalism and results presented 
in this paper are unprecedented in this respect and are crucial 
steps toward applying the BAO measurements in redshift sur- 
veys to constraining the nature of dark energy, the curvature 
of the Universe, etc. 

Considering the biasing effects in galaxy clustering is be- 
yond the scope of this paper. A simple prescription of the 
biasing is to replace f ^ /3 - f/b in our equations, where b 
is the linear bias parameter However, nonlinearity and nonlo- 
cahty of biasing might not be negligible even on BAO scales 
lfTolll2[l . Incorporating the nonlinear and nonlocal biasing into 
our formalism will be addressed in future work. 
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FIG. 9: Nonlinear evolution of the baryon acoustic oscillations in redshift space for various redshifts, z = (top left), 0.5 (center left), 1 
(bottom left), 2 (top right), 3 (center right), 4 (bottom right). Direction dependencies are spherically averaged. Each power spectrum is divided 
by a smoothed, no-wiggle linear power spectrum P„„(k) in real space |44]. Black (solid) line: this work; red (dotted) line: linear theory (Kaiser 
formula); green (short dashed) line: 1-loop SPT; blue (long dashed) line: a simple model of redshift-space power spectrum t35i1 . The validity 
range k < ^nl/S, where our result is expected to be accurate within a few percent, is indicated by arrows. 



APPENDIX A: INTEGRATING CUMULANTS OF THE 
DISPLACEMENT FIELD 



In this appendix, we present results of integrating cumu- 
lants of the displacement field, which appear in Eq. ( |35] ). Al- 
though we do not need to evaluate the integrations as tensors 
in real space, we do need to do that in redshift space. 

To present the results, it is convenient to define the follow- 
ing integrals: 



where 



Qn(k) 



R„(k) 



1^3 noo 



xPlWI -2rx)^^^] 



Qnir, x) 



1 k^ 



(1 + -2rx)2' 
drPakr)R„(r), 



(Al) 
(A2) 



e, =r2(l-x2)2, 

22 = (1 -x^)rx(l - rx), 

Q3^x\\-rx)\ 

24 = 1 - x\ 



and 



1 +r 



/?2 = ^(l-r')(3-2r^ + 30 



l + r 
l-r 
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FIG. 10: Comparison of the power spectra to the A'-body simula- 
tions of Ref. 1 8] in redshift space at redshifts z = 0.3, 1, and 3 
from top to bottom, respectively. Open circles: A'-body simulations; 
Black (solid) line: this work; red (dotted) line: linear theory; green 
(dashed) line: 1-loop SPT; blue (long dashed) line: a simple model of 
redshift-space power spectrum 13511 . Only nonlinear deviations from 
the linear growth are measured in A'-body simulations to reduce cos- 
mic variances. Estimates of validity range, k < are indicated 
by arrows. In the z = 3 sample, the range where contributions from 
the naive shot noise exceed 1% is also indicated. 



The following expressions are directly obtained by substi- 
tuting Eqs. (|23ll-(|26ll into Eqs. (EB-®: 









(A9) 


C<f (fc) = 


98 ^ 




(AlO) 




u 21 




(All) 



because the rotational covariance implies 



I 



(In} 



■g{k, p)T{k, -p, p) oc k. 



(A12) 



where g is a scalar function. 

The remaining integrals are obtained by using the rotational 
covariance. For example, when the integrand is given by pipj 
times a scalar function of fe and p, integration over p should 
result in a form XSij + Yk/kj. The coefficients X and Y can be 
obtained by contracting the original expression with 5,^ and 
kikj, and so forth. After these kinds of manipulations, we ob- 
tain 



d^P ^(112)., , X 

C;,,, (k,-p,p- k) 



(2;r)3 ""'V* 



3 kik jkf^ 



&}^^\k,-p,p-k) 



14 k^ 



3 ki6jk 



(Pp 



Cf'\k,-p,p-k) 



(2;r)3 

3 kik jkif 
~ 14 k(> 



(Qi + 2Q2) ■ 



14 k^ 



3 kjSjk 
14 k* 



Qi 



(A13) 



(A14) 



r ^c;'.'>(p)d">(fe-p) 



224) 



1 kiikikkkt) 

+ g ' \ 3Qi+ 2462 + 863 -464), (A15) 



In deriving Eq. ( lAl lb . the transverse part of Eq. ( |25] | vanishes. 



where the spatial indices are symmetrized over in Eq. ( IA15I) . 

If we substitute those expressions of Eqs. ( |A9I )-( IATT] ) and 
(IA13l l- (IA15b into Eq. ( |35] ), the real-space result of Eq. ( |36] | is 
reproduced. 



APPENDIX B: SPT CALCULATION OF THE 1-LOOP 
POWER SPECTRUM IN REDSHIFT SPACE 

In this appendix, the nonlinear power spectrum in redshift 
space is calculated b y 1- loop SPT. A formal expression found 
in the literature l(33l l34ll is given by 3-dimensional integrals, 
which are reduced to lower dimensional integrals below. It is 
suggested that the 1-loop calculation of the power spectrum 
in redshift space by SPT does not give satisfactory results, be- 
cause SPT in redshift space breaks down at larger scales than 
in real space |34]. In this paper, we do not recommend the 
extensive use of the 1-loop SPT in redshift space. However, 
it is still useful to calculate the nonlinear corrections by SPT 
in redshift space for the following reasons: First, the 1-loop 
corrections indicate when the linear theory is applicable and 
when nonlinear effects become important in redshift space. 
Second, the 1-loop SPT gives the asymptotic behavior of non- 
linear corrections in the ^ — > limit, which should be com- 
pared with any nonlinear theory of redshift-space distortions 
on large scales. 
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FIG. 11: Nonlinear evolution of the baryon acoustic peak in redsiiift space for various redshifts, z = (top left), 0.5 (top rigiit), 1 (bottom left), 
2 (bottom right). Angular dependencies are averaged. Black (solid) line: this work; red (dotted) line: linear theory (Kaiser-Hamilton formula). 
green (dashed) line: Gaussian-filtered linear theory (see text). 
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FIG. 12: Comparison of the correlation functions to the A'-body sim- 
ulations of Refs. (si l46ll in redshift space. Open circles: A'-body 
results; Black (solid) line: this work; red (dotted) line: linear the- 
ory. Only nonlinear deviations from the linear growth are measured 
in A'-body simulations to reduce finiteness effects. 



The comoving real-space position x and the comoving 
redshift-space position s are related by Eq. ( l44l l. or 

s - x + z. (Bl) 

aH 



The density field in real space p(x) and that in redshift space 
Pj(s) are related by conservation relation: 

Ps(s)d^s = p{x)cPx, (B2) 

Therefore, the density contrast in redshift space is given by 

5s(s) = [1 +^(a;)]7"' - 1, (B3) 

where J - d{s)/d(x) is the Jacobian of the mapping from real 
space to redshift space. One can easily calculate the Fourier 
transform of Eq. ( IB3l i: 

6,(k)^6(k) + J d\e-'''''(e'''-'"^ -l][l+6(x)], (B4) 

where u~ = -v~/aH. This relation is applicable even in the 
fully nonlinear regime. In Fourier space, M,(fe) - ikfi{k)/k^. 

Following the route of SPT, the density contrast and the pe- 
culiar velocity, which are absent in an homogeneous universe, 
are expanded in perturbative series: 

e = 6»*'* + 6i*2> -I- 0^^' + ■■■ , (B6) 

where = V ■ v/aH corresponds to the velocity divergence. 
In the SPT, only growing-mode solutions in each order are re- 
tained, and the peculiar velocity field is consistently assumed 
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to be irrotational. Thus, the velocity field is fully character- 
ized by the velocity divergence |[20.1 . In Fourier space, each 
perturbative term is given by 



P± 



X Gnipi 



...,p„)6o(P]) 



k 



V 1=1 



■5o(p„), (B7) 



(2nY 

■ ■,P„)5q{pi)- 



Vi=I 

SaiPn). 



(B8) 



in real space ifTsifTgl] . The power spectrum P^ik) in redshift 
space is defined by 



(6,{k)Uk'))^ = {2nf6\k + k')P,{k). (B13) 



Note that the power spectrum in redshift space is no longer 
isotropic. The perturbative expansion of the power spectrum 
up to second order in the linear power spectrum (i.e., fourth 
order in (Jo) is given by 



The perturbative kernels F„ and G„ do not depend on time in 
an Einstein-de Sitter model. They weakly depend on time and 
cosmological parameters in general cosmology, but it is still a 
good approximation that the perturbative kernels are replaced 
by those in the Einstein-de Sitter model even in arbitrary cos- 
mology lESl . The expUcit form of the perturbation kernels is 
given in Refs. lliiolill- 

Expanding the peculiar velocity field in Eq. ( IB4b . and ap- 
plying the perturbative expansions of Eqs. ( IB5I ) and (IB6I) . we 
obtain the perturbative series for the redshift-space density 
field: 



X S„{pu. 



j1 In 
,Pn)(5{)(pl)-- 



(B9) 



where the perturbative kernels in redshift space S „ are given 
by 



PAk) = D\t)P,n(k) + D\t) [Psiiik) + P,n{k)] , (B14) 



where P^n is the contribution from P^22 is from 

{WW^), and f si3 is from + The first quan- 

tity is the linear (or tree-level) power spectrum, and the last 
two quantities are 1-loop corrections, which are given by 



P,„(fc)= [5i(fc)]V„(^) 

V3 



(B15) 



P.ii{k)^2 j ^[si{p,k-p)fPo{p)Poi\k-p\), 



(B16) 

P,i3(fe) = 65,(fc)Po(fe) r ^S\(k,p,-p)PQ(p), (B17) 
J (2ny 



5i(pi) = 1+//, 

SliPuPl) = F2(pi,P2) + fl^^G2(Pl,P2) 



(BIO) 



2 2 
Pi P2) 



^ 1.., .2P}zP2z 

2 P\pI 



(BID 



Si{p\,P2,Pi) = Fi(puP2,P3) + ff^ G3(pi,P2,Pi) 



+ fkfI^F2(p2,Pi) + fkjjp^^—^G2(P2,P3) 
Pi 



T 

Ip2 + P3r 



+ (fkfi) 



: Plz(P2z + P3;) 
p1IP2 +P3P 



G2(P2,Pi) 



^ ^.fl ^2 PlzP2z 

2 P\pI 



1,., .3PlzP2zP3z 
Tifkp) 2 2 2 ' 

6 pApI 



(B12) 



up to third order, where ^ - k ■ z/k, k = pi + ■ ■ ■ + p„, and 
pi^ - Pi ■ z etc. The expression of 5 3 should be symmetrized 
over its arguments when necessary. These kernels are equiv- 
alent to those derived in Ref. [34], while the derivation and 
apparent expressions here are somewhat different from those 
of Scoccimarro, Couchman, and Frieman. When the linear 
bias b is present, the perturbative kernels for density fluctua- 
tions of galaxies are given by replacing f ^ fi - fjb m the 
above kernels. When the bias is modeled by a nonlinear and 
local one, it is straightforward to include that by a method of 
Ref. m. 

Given the perturbative kernels, the nonlinear power spec- 
trum in redshift space is calculated in a similar way to that 



where 5 is a symmetrized kernel obtained from S „ . The first- 
order power spectrum is given by 



Psii(fc) = (l +fp^fPn{k) 



(B18) 



which is the Kaiser formula H31I1 in linear perturbation theory. 

Integration over the azimuthal angle can be analytically per- 
formed in Eqs. ( IB 16b and ( IB17I ). and further integration over 
the polar angle in Eq. ( IB17b is also possible. Terms which 
contain the line-of-sight components of the wave vector can 
be evaluated by using the rotational covariance as described 
in Appendix lAl After a lengthy but straightforward calcula- 
tion, the final results are given by 



Ps22(k)^Yj^'^"f"'^ drPoikr) j ^dx 



{1+ r^- 2rx) 



2' 



(B19) 



P,13(fc) = (l+/A^Vo(^) 



YJ^'''"f'"^ Jo drPo(kr)B„„,(r) (B20) 
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where non-vanishing components of A„,„ and B„,„ are 



and 
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—(-7 + 35x^ + 54-rx - 1 lOrx" 
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A34 = -(1 - x^)(2 - 3r^ - Urx + 15rV), 
A44 



—(-4 + 12x2 ^ 3^2 ^ 24rx 
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(B21) 
(B22) 
(B23) 



+6r2(19 - 157x2 + 236x4)] , (B24) 

(B25) 
(B26) 



(B27) 
(B28) 



- 30r^x^ - 40rx^ + 35r2x'*), (B29) 
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(B30) 
(B31) 

(B32) 

(B33) 
(B34) 



In a limit / — > 0, only m - n = Q terms survive, and the 
previous results of the 1-loop power spectrum in real space 
il8ll are recovered. 
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